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^ ' Recent experimental measurements and lattice QCD calculations are now reaching the 

^^ , precision (and accuracy) needed to over-constrain the CKM parameters p and f). In this 

l/^ ■ brief review, I discuss the current status of lattice QCD calculations needed to connect the 

^S) ' experimental measurements of B meson properties to quark flavor-changing parameters. 

^+ ' Special attention is given to _B — > rriv, which is becoming a competitive way to determine 

C^ , |Vut.|, and to B^ — B^ mixings, which now include reliable extrapolation to the physical 

VP 1 light quark mass. The combination of the recent measurement of the Bs mass difference 

f**^ ' and current lattice calculations dramatically reduces the uncertainty in \Vt^\. I present 

an outlook for reducing dominant lattice QCD uncertainties entering CKM fits, and I 

remark on lattice calculations for other decay channels. 
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1. Introduction 



, . , The main motivation for studying the properties of B mesons is to understand 

d ' how quarks change their flavor. Is the Cabibbo-Kobayashi-Maskawa (CKM) mech- 

anism enough to describe flavor-changing processes, or can we see the effects of new 
physics? 

Even with the extraordinary effort made by the B factories, we presently see no 
cracks in the CKM model. This allows us to eliminate or tightly constrain many new 
physics models, afflicting them with a so-called flavor problem. Continued agreement 
of experiment with the CKM model will force us to explain the minimal flavor 
violation. Better yet, fissures will appear and possibly favor one new model over 
another. As usual in particle physics, this type of indirect search complements and 
guides direct searches. 

Since quarks are confined inside hadrons, precise and accurate theoretical calcu- 
lations are necessary to connect the experimental measurements of meson decay and 
mixing to the fundamental quark couplings. Many quantities important for studying 
quark flavor are calculated using lattice QCD (LQCD). This brief review focuses on 
calculations of B decays and mixings. Sec.|21is an overview of the methodology. 

Semileptonic decays like B — > irlv and B — > Div provide information about 
how the h quark decays to a u and c quark by W emission. Since these decays 



February 2, 2008 5:52 WSPC/INSTRUCTION FILE wingate-mpla 



2 Matthew Wingate 

proceed doniinantly through tree diagrams, they aUow one to determine the CKM 
matrix elements \Vub\ and \Vcb\, given clean experimental measurement of decay 
rates and reliable theoretical calculation of hadronic matrix elements. B — > niu 
decays were hard to use to determine \Vub\ because of small branching fractions 
and large theoretical uncertainties. In the past year both experiment and theory 
have improved greatly, the latter through unquenched LQCD calculations. Sectional 
reviews this important development, gives an update on _B — > Dliy, and discusses 
the difficulties and possibilities for lattice calculations relevant to other semileptonic 
decays which might also determine | Vub \ ■ 

Mixing between neutral mesons proceeds through loop diagrams to which new 
physics particles can contribute. It is possible that new physics could be discovered 
by assuming that only Standard Model particles contribute and showing that the 
resulting CKM parameters do not agree with the CKM parameters obtained from 
tree-level decays. Both K^ — K^ and i?° — B^ oscillations provide important ways 
to search for deviations from the CKM mechanism. As in the case of semileptonic 
decays, calculations of hadronic matrix elements are necessary in order to learn 
about quark level processes from experimental measurements of mass and lifetime 
differences. Section 0] describes recent progress made using LQCD in concert with 
new experimental results. 

Rare decays, like B -^ if*7, occur through loop diagrams called penguins. In 
principle lattice QCD can also provide reliable calculations of the relevant hadronic 
matrix elements for these processes. Some technical difficulties facing LQCD con- 
tributions to rare decays are explained in Section [5| 

This review closes with the lattice prediction and subsequent experimental mea- 
surement of the Be mass (SecEJ and a few concluding remarks (Sec.|7|). 

2. Overview of Methodology 

To many people lattice QCD is a black box into which one puts the strong coupling 
constant and quark masses and out of which one obtains the properties of hadrons. 
If only that were so. Those people who live inside the box see LQCD is a method for 
numerically integrating Euclidean spacetime path integrals, valid even when there 
are no small parameters about which to form a perturbation theory. Communication 
between those living inside and outside the box is imperative if progress is to be 
made. Here I discuss some details of LQCD that will facilitate later discussions 
about recent progress and future improvements and obstacles. Interested readers 
can find more detailed treatments elsewhere.^ 

Expectation values of observables are given by integrals over all possible field 
configurations. For example the expectation value of operator O is given by 

(O) = 1 /"l?V^V^^C/0[V',^,t/] e-^I'^-'^'^] (1) 

where the tp, 'tp represent the quark, antiquark fields and U represents the glue field. 
The action S is the integral over space and imaginary time of the QCD Lagrangian. 
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Lattice QCD provides a method for numerically evaluating path integrals from first 
principles. The integral is made finite by discretizing spacetime. The presence of 
the damped exponential in |^ means that only a small set of field configurations 
will contribute to the path integral. This is analogous to studying a classical gas: 
the important configurations of particle positions and velocities are a small fraction 
of the multitude of possible configurations, most of which give exponentially small 
contributions to statistical traces over microstates. 

The way that fermion antisymmetry is represented in path integrals causes great 
difficulty for numerical methods. The quark and antiquark fields are valued over the 
field of anticommuting complex numbers (Grassmann numbers) instead of the usual 
complex numbers. This is a formal construction which generally does not permit 
numerical evaluation. Luckily, if the action is quadratic in the fermion fields, Sf = 
ipQ^p, the integral over fermion fields can be done exactly, yielding a determinant, 
detQ. For example, the vacuum expectation value of a quark bilinear is given by 



(iPVi;) = /p[/ ArAe-fQ"'[^lC detg[[/]e- 
J oC oC 



■s,iu] 



(2) 



where F is a product of Dirac 7 matrices and the C, C act as fermion sources. The 
quantities calculated for B physics results, 2- and 3-point correlation functions, are 
straightforward generalizations of (0) • The penalty we pay for this formal integration 
is that a nonlocal updating algorithm must be used to do importance sampling. A 
necessary step involves inverting a matrix which becomes singular as up/down quark 
masses are decreased from artifically heavy to their physical values. 

Writing ^ in detail allows us to define some terminology. Note the 2 places 
the quark matrix Q appears in ^ . These correspond to 2 distinct steps in a lattice 
QCD calculation. The determinant of Q is included in the probability weight during 
the generation of typical, important glue field configurations. The resulting collec- 
tion of configurations is then stored and later used to calculate ensemble averages 
of operators. In the case of fermionic operators, the inverse of Q is computed during 
this ensemble averaging. The term "sea quarks" refers to contributions of det Q to 
the importance sampling and "valence quarks" refers to the fully dressed propaga- 
tors Q~^. Logistically nothing prevents one from using different quark mass values 
for the sea and valence quark steps. This procedure is called "partial quenching" 
and is useful since the generatio n of configurations is much more costly than the 
calculation of quark propagators.'^'"' Full QCD is in some sense a special case of 
partially quenched QCD where irisea. = "^valence- 

The last important point is that the discretization Q of the Dirac operator, 
Ip + m, is not unique. The choice of discretization has consequences for the sys- 
tematic uncertainties in numerical results and for the computational expense. One 
discretization, the improved staggered quark action, permits numerical calculation 
with sea (and valence) quark masses down to 1/10 the strange quark mass. (The 
cost of the calculations increases rapidly as m decreases.) Ensembles of configura- 
tions which include effects of 2 + 1 fiavors of sea quarks using an improved staggered 
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action, have been generated and made public by the MILC collaboration!^ The ef- 
fects of using light sea quarks are dramatic: cleanly computable quantities, many 
of which the quenched approximation got wrong by over 10%, agree nicely with 
experiment!^ Calculations with other quark discretizations are desirable (1) to fur- 
ther check the "fourth-root" algorithmic trick one needs with staggered quarks to 
obtain 2-1-1 sea quark flavorsPl (2) to avoid the complications staggering creates in 
calculations of light baryon properties.'-'J^'^Unquenched calculations with improved 
Wilson quarks and domain wall quarks are making steady progress. However, their 
advantages are tempered by having to extrapolate using data with heavier quark 
masses, perhaps without sufficient theoretical control. Checks similar to those done 
for staggered fermions are needed for non-staggered actions to demonstrate that the 
extrapolations to physical sea quark masses are under control. 

3. Semileptonic Decays and CKM Matrix Elements 

Semileptonic decays, where quark flavor is changed by emission of a single W boson, 
provide direct access to the first 2 rows of the CKM matrix. It is useful to focus on 
these tree-level decays separately from processes which begin at 1-loop level because 
new physics particles generally should affect the latter more than the former. 

In order to use experimental measurements of semileptonic decays to determine 
the CKM matrix elements, one needs accurate calculations of the weak current be- 
tween hadronic initial and final states. In the next section the lattice calculations 
needed for \Vcb\ and \Vub\ are discussed (as are \Vcs\ and \Vcd\, briefly). The mixing 
calculations relevant for \Vtd\ and \Vts\ are discussed in Sec.0I This brief review will 
not cover calculations for \Vus\ using lattice QCD, although there has been signif- 
icant recent progress, using the leptonic decay constants " --^J or the semileptonic 
form factor! I I I I Finally, \Vud\ is most preci sely determined from nuclear and 
neutron beta decays without needing lattice QCD.l^ 

3.1. B -^ Z)(*)£i/ 

Presently, the most precise determination of \Vcb\ comes from inclusive 
measurements.!^ Improvements in lattice QCD calculations will help to reduce 
the uncertainties in determinations from exclusive decays, especially B -^ D*(.v. 
Regardless of which method is most precise, cross-checks are important to truly 
over-constrain the CKM model. 

The matrix element \Vcb\ may be determined from exclusive semileptonic B 
decays to the D or D* by fitting the experimental measurements of the differential 
partial decay width to 

^[B^ Di*Hu) = ^|i^/C^^^<., [w)Ti^^,„ {w) (3) 

where w ^ vb ■ Vjj(*) is the velocity transfer and /Cs^^c*) [w) is a known kinematic 
function of w. The shape of J-'b~>d(*) i'w) near w = 1 is given by dispersive bounds 
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and HQETJ^and experimenters quote measured values for Tg^ij(,){l)\Vcb\- The 
overall normalization can then be computed using lattice QCD. 

The LQCD calculation can be made very precise by constructing a double ratio 
of matrix elements, wherein many uncertainties cancelf**' Preliminary unquenched 
results give JFb^d(I) = 1.074(18)(16), where the first error are statistical and the 
second is a combination of 1% systematic errors M"! Combining this with the Winter 
2006 experimental averagJ22 gives \Vcb\ = (39.7 ± 4.2|exp ± 0.9|iatt) x 10"^. 

While for B —^ Dip the dominant uncertainty lies on the experimental side, 
the experimental error is 2% for B -^ D*£v. Furthermore, Luke's theorem says the 
heavy quark effective theory errors should also be smaller for the latter decay. One 
hiccup is that the LQCD extrapolation oi Tb^d* (1) from the input light quark mass 
to the physical light quark mass passes through the D* -^ Dir threshold. Improved 
staggered fermion calculations are poised to solve this problem for 2 reasons. First, 
the calculations are able to be performed within the chiral regime, so that heavy 
meson chiral perturbation theory can provide a reliable formula for extrapolation. 
Second, the presence of artificially heavy staggered pions actually smoothes out the 
cusp in J^B^D'i^) at threshold as seen in staggered chiral perturbation theory!^ 

3.2. B -^ nlv 

Because tt and B have the same parity, only the vector part of the weak current 
mediates B — > ttIp. The matrix element can be parameterized by 2 form factors: 

(nip^WnBipB)) = .Mq') (p'b+p'. "'^^~"^ ^) + /o(g^) "'^^~ "' g" • (4) 

The form factor f+{q^) describes the momentum dependence of the differential 
partial decay rate in the B rest frame by 

^ - ||3KnK.n/+(.^)p. (5) 

Experiments have measured the B^ -^ n^l^iy branching fraction with impressive 
precision. \Vub\ is determined by integrating /+ over q^. There is presently a lower 
limit on q^ below which the lattice calculations incur large discretization errors: 
the TT should not have momentum larger than the inverse lattice spacing. Therefore, 
theoretical errors are minimized by integrating over 1 6 Ge V < q^ < (/max and using 
the branching fraction with the same momentum cut.'^"' 

Two groups are completing calculations of the B — > -kHu form factors. They 
both use the MILC configurations, which include effects of 2+1 quark flavors using 
the improved staggere d ac tion. The main difference is the heavy quark formulation: 
either lattice NRQCeB3I or the Fermilab method. -'^^ NRQCD is computationally 
inexpensive but requires the heavy quark mass to be greater than the inverse lattice 
spacing, precluding charm studies with NRQCD on the MILC lattices. The Fermilab 
method is valid for all quark masses. 
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Fig. 1. HPQCD results for the B ^ tt form factors. The symbols are lattice data after extrapo- 
lating the light quark mass to zero, and the curves are fits to the Ball-Zwicky parameterization. 



One needs to fit the q^ dependence of the form factors to a functional form 
for two reasons. First, to interpolate the finite quark mass lattice data to fixed 
values of i?^, so that the chiral extrapolation may be performedM^I Second, to 
facilitate integration of /-f over a range of q" corresponding to the same range for a 
measured branching ratio, yielding a det ermination of |y„b|. The commonly- used 3- 
parameter Becirevic-Kaidalov ansata^fits the HPQCD lattice data well only after 
extrapolation to the chiral limit. Instead we use the 4-parameter Becirevic-Kaidalov 
ansatz advocated by Ball and Zwicky.!^ 

The final HPQCD result iJ 



dq- 



dV 



1.46 ± 0.23 ± 0.27 ps' 



(6) 



iKbP Jie Gov^ dq^ 
where the first error is statistical plus chiral extrapolation and the second is 
other systematics added in quadrature. The dominant source of systematic un- 
certainty is due to truncating at 1-loop order the perturbative matching between 
the lattice NRQCD current and the continuum weak vector current. This agrees 
within errors with the preliminary result using the Fermilab heavy quark action: 
iKfor^/ifc^v^ dq^{dT/dq^) = 1.83 ± 0.50 ps-^T^Thc dominant error in the lat- 
ter calculation is due to discretization errors in their heavy quark action. Taking 
curren t ex perimental results for the i?° — > 7r~£"'"^ branching ratio °" - ■■^ and B^ 
lifetimtEl along with the HPQCD result ® yieldJ^ 



\Vub\ = (4.22 ± 0.30exp ± 0.34,tat ± 0.38sys 



) X 10" 



(7) 



With NRQCD for the heavy quarks, the leading uncertainty comes from the per- 
turbative matching of the lattice currents to the continuum weak current!^ The 
matching is presently done at 1-loop level. The 2-loop calculation is very involved 
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and will require automated perturbation theory techniquesM^ With the Fermilab 
heavy quark formulation, there is a conserved current which allows most of the 
matching to be determined nonperturbatively. The remaining perturbative correc- 
tions are tiny, so truncating at 1-loop level is sufficient until other errors are de- 
creased. The leading uncertainty in this case is the discretization error in the heavy 
quark formulation. The improved action has been worked out, but an efficient 
implementation is still under development. 

There are several efforts to reduce any systematic uncertainties due to model- 
ing the shape of the form factor. The shap e is str ongly constrained by unitarity, 
analyticity and heavy quark effective theory! '^ 1 ^*^ I Plotted against an approp riate 
variable, the shape of /+ is consistent with a linear fit .1^ Recent analyse d I I make 
it apparent that LQCD calculations can reduce the uncertainty in \Vub\ by more 
precisely determining the form factor normalization. 

We see that determinations of \Vub\ from B -^ tt£i/ decays are becoming com- 
petetive with the precision from inclusive determinations.!^ Using exclusive decays 
also has the advantage of avoiding the hypothesis of quark-hadron duality which is 
necessary for inclusive determinations and blunts tests of the CKM model. 



3.3. D -^ Tziu and D -^ K£u 

Even though D decays are outside the scope of this brief review, the recent calcula- 
tions of semileptonic D form factors have implications for the B ^> tt form factors, 
on top of providing determinations of \Vcd\ and | Ks | ■'"'^ The lattice calculations are 
done with a subset of the same MILC configurations discussed above^ and use the 
improved staggered quark action in both the B and D correlation functions!"^ The 
charm quark is treated using the Fermilab formulation, as in the Fermilab/MILC 
B meson calculations. Given these similarities, the agreement between the l attice 
calculation of f+{q^) for D ^ K and the recent experimental measuremenl! ^^ | 4i | 
testifies to the validity of the B ^ tt form factor calculations. 



3.4. Other b — > uiv Decays and Associated Challenges 

Exclusive determinations of \Vub\ from other exclusive B decays B -^ t:(.v are also 
valuable. Every independent measurement has the potential to play a role in the 
quest to see new physics through over-constraining the CKM parameters. Even 
within the Standard Model, it is conceivable that another decay channel will lead 
to a more precise determination than B -^ nlv, provided that lattice QCD can 
meet some challenges. 

Experiments show the B —^ p(!.v branching fraction to be twice as big as to -niv^ 
so its precision is easier to reduce. However, to get an accurate \Vub\ the theoreti cal 
uncertainties in the form factors must be under control; presently they are not!^ 
The main open problem is the exptrapolation of lattice data obtained with quark 
masses where the p is stable through p -^ tttt threshold. The firm theoretical ground 
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provided by chiral perturbation theory for B and n properties becomes quite muddy 
underneath the p. ^^ 

Experiments have measured B — > rjiiy}^ The lattice calculation with the 77 
suffers because it is partially a flavor singlet. There are disconnected contributions 
to the correlation functions where the valence qq pair are annihilated then recreated. 
These are notoriously hard to calculate due to poor signal-to-noise. Nevertheless, 
these contributions are important to the 77 propagator du e to the anomaly. 

The situation may not be as glum for B — > usli/^^ which has recently been 
observed!^ The w is a narrow resonance in contrast to the p. Despite being a 
flavor singlet, violations of the OZI rule are expected to be smaller for the vector 
w than for the pseudoscalar 77, since the anomaly affects only the latt er. LQ CD 
calculations of flavor singlet meson masses support this expectation so far! ^ I ^i l'^y^ 
should explore in more depth the reliability of calculating the properties of the lo 
on the lattice and see how far we can push a calculation of the B ^f u) form factors. 
If we can truly neglect disconnected contributions to the uj propagator, then the 
lattice calculation of B ^ a; form factors would be identical to B ^ p, with less 
headache regarding the meson's width. Disconnected contributions are routinely 
ignored in LQCD calculations of the (j) meson mass, for example. 



4. Mixing and the Search for New Physics 

Studying the oscillations of neutral mesons has the potential to expose physics be- 
yond the Standard Model in a way that the semileptonic decays do not. Any new 
physics entering the semileptonic decays would be a small contribution to the tree- 
level standard model decays. In contrast oscillations begin at one-loop level in the 
standard electroweak theory, and new particles can enter directly in those loops. 
The leading uncertanties so far are theoretical. Nevertheless, with increased preci- 
sion in lattice calculations and with the new measurement of B'^ — B^ oscillations, 
the potential to see a breakdown of the CKM model is realistic, even if unreal- 
ized presently. This section summarizes the most recent lattice QCD results and 
prospects for further improvement. 

Before turning to the B" mixings, let us consider K'^ and D'^ mixings. Neutral 
kaon mixing is another important window into flavor physics which relies on lattice 
QCD. The theoretical error presently dominates that constraint on (p, fy). This 
situation should improve as we see mor e unquenched calculations of the hadronic 
matrix element parametrized by Bk^^ 

£)0 _ jjo jnixing proceeds very slowly compared to mixing the K and B systems. 
In the Standard Model, this is because the top quark mass, more than 40 times mi,, 
enhances the AB — 2 and AS* — 2 box diagrams compared to AC — 2. Of course, 
rare processes are just the kind one wants to study to discover new physics. As 
soon as precise experimental measurements can be made, lattice QCD would be 
valuable in connecting the meson mass and width differences to models of quark 
flavor-changing interactions. The relevant matrix elements were calculated some 
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Fig. 2. Ratio of /a^ //g {times , I mg^/rriB ) vs. mass of the light quark (in units of the physical 
strange quark mass). 

time ago '^ l ' ' *'^ ! and should be updated using unquenched configurations and the 
Fermilab formulation for the charm quark. 

4.1. B° -W 

In the Standard Model 5° — B^ mixing proceeds via 1-loop box diagrams, with the 
dominant contribution coming from a top quark in the loop. Measurements of the 
oscillation frequency Anid can be used to constrain the CKM matrix element \Vtd\- 
The relevant hadronic matrix element is conveniently parameterized as 



{B°\ibd)v-Aibd)v-A\B'^) = g/ImiBB . 



(8) 



This is helpful for two reasons. The leptonic decay constant /b is simpler to com- 
pute in lattice QCD. In addition chiral perturbation theory shows that fs is more 
sensitive to the quark mass than Bs, the ratio of the full matrix element to its value 
in the vacuum saturation approximation. 

As discussed above, the LQCD calculations are done with quark masses mq 
larger than the physical up/down quark mass. Let's denote the corresponding decay 
constant by /b ■ In the past, the extrapolation of /s to the physical light quark 
mass was done phenomenologically. Since the lattice data for Jb are linear in 
the regime ms/2 < iriq < 2ms, a linear fit was used for the extrapolation. Chiral 
perturbation theory gives the ruq dependence of /b including logrithmic curvature; 
however the theory must be used in the small mq regime, approximately mq < ms/2. 
The use of improved staggered quarks for the fermion determinant and for the 
light quark propagator allows LQCD calculations to be performed inside this chiral 
regime. 
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The HPQCD collaboration has used the MILC lattices discussed above to cal- 
culate Jb vs. mq (see Fig. |21). The data are fit using an extension of chiral per- 
turbation theory suitable for staggered fermion calculations.!^ Our result for Jb is 
216±9stat±20sys MeV.I^The leading uncertainty is due to the perturbative match- 
ing between the currents in the heavy quark effective theory (lattice NRQCD) and 
the continuum weak axial vector current, which has been done at 1-loop order. 12^ As 
in the case of the _B — > tt form factors, reducing this uncertainty requires a difficult 
2-loop calculation. 

Very recently the Belle experiment reported a result /s = 176(t23)stat(tig)sys 
MeV from a first measurement of B^ -^ T~Vr combined with |Vu(,| from B -^ 
Xulv. ^\\% confirmation of the LQCD prediction directly speaks to the validity 

of LQCD calculations for B meson properties, in particular the use of NRQCD heavy 
quarks, staggered light quarks, and staggered chiral perturbation theory. 

In order to further judge the reliability of the lattice QCD calculation of /s, 
it is useful to look at the D decay constant. A lattice calculation, which differed 
from the /s calculation only in heavy quark action, recently predicted /n to be 
201±3stat±17sys MeV,'^much more precise than the 20% experimental uncertainty 
at the time.ESl Subsequently, CLEO-c found f^^^ = 223 ± IT^tat ± %y^ MeV from 
£)+ — > \i^v'^^ This successful prediction by lattice QCD using the same MILC 
configurations, light quark action, and staggered chiral perturbation theory provides 
a significant test of the Jb calculation. 

The remainder of the hadronic matrix element, B^, remains to be calculated on 
the MILC lattices, although preliminary efforts have begun. '^ The state-of-the-art 
is still the JLQCD calculation done with improved Wilson fermions.''^ The same 
systematic uncertainties affect their result as the HPQCD result for Jb'- truncating 
the perturbative and heavy quark expansions. Since the chiral extrapolation for B^ 
is milder than for /s, it is reasonable to expect their result to be close to what will 
be obtained with ligh ter quark mass computations. Combining the 2 groups' results 
one obtains /sv^B^ = 244 ±26 MeV.I^This reduces the uncertainty down to 11% 
from the previous world average JbV^b = 214 ± 38 MeV.I^ 

Fig. Ot shows the constraint in the p — fj plane from Atti^. The improvement 
indicated by the 2 bands is achieved by reducing the long extrapolation of /b in 
light quark mass. 

For physics beyond the standard Model, B^ — B'^ oscillations may be governed by 
4-quark operators with different structure than in JHI); this is certainly the case for 
supersymmetric extensions of the Standard Model." Therefore, matrix elements of 
all 5 possible operators should be computed using unquenched lattice QCD, as has 
been done in the quenched approximation.l^^ 

Given the difficulty of 2-loop calculations in lattice perturbation theory, 
compounded by complexity of heavy quark lattice actions, a fully nonpertur- 
bative met hod would be most welcome. Several such approaches are being 
investigated! I I I I Only time will tell which of these approaches will first sur- 



February 2, 2008 5:52 WSPC/INSTRUCTION FILE wingate-mpla 



B PHYSICS ON THE LATTICE: PRESENT AND FUTURE 11 



(a) (p, r|) from Am^ only 

T — ' — I — ' — I — ' — I — ' — I — ' — I — ' — r 



(b) (p, r|) from Am^ and Am^ 





Fig. 3. Constraints (at the 2a level) on \Vtd\ = A\^^{1 — p)^ + fj"^. (a) Old and new constraints 
from Arra^j. (b) Constraint from Am^j/Ams using the new CDF measurement and recent LQCD 
improvements. 



pass in precision the methods used in present calculations and when that might 
occur. 



4.2. B° - BO 

s s 

Unlike leptonic B, D, and Ds decays, which are merely helicity-suppressed, lep- 
tonic Bg decays are forbidden in the Standard Model since they would pro- 
ceed through a flavor-changing neutral current. Nevertheless, the decay constant 
/b^ = {0\AQ\Bs)/mB^ is finite and approximately equal to /b- Again, it provides a 
useful parameterization of the hadronic matrix element relevant for B^ — B^ mixing, 
as in ||SJ). 

Very recently, the DO experiment reported a 2-sided bound on t he oscillation 
frequency 17 ps~^ < Amj < 21 ps^^ at the 90% confidence levelp^and the CDF 
collaboration reported a measurement of Ato^ = 17.3 f^'\ ps^^. " Some theoretical 
uncertainties cancel in ratio of the neutral Bg and B oscillation frequencies, so a 
theoretical calculation of Jb ^bJ Ib^b can now provide a much tighter constraint 
on \Vtd\ than using only nonstrange B oscillations. 

The HPQCD result for /^^ is 260 ± 7stat ± 28sv s MeV using a nonrelativistic 
b quark and the MILC lattices described above .1^ As with Jb, the perturbative 
matching is the dominant source of systematic uncertainty. This uncertainty largely 
cancels in the ratio JbJIb = 1-20 ± 0.03 ± 0.01, where the first uncertainty is 
statistical error in the correlation functions and in the chiral extrapolation, and the 
second uncertainty estimates truncated terms in the usual expansions!^ Combining 
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this ratio with the JLQCD ratio of B factors^Q yieldJ^ 




1.210 1^^ (9) 



Ratios of the B and Bg decay constants and B factors have also been determined in 
the static hmit with 2 flavors of sea quarks using the domain wall fermion action!^ 

The constraint obtained by combining the improved LQCD ratio l|njl with the 
new CDF measurement"*^ (and A = 0.225(1) ) more than halves the width of the 
allowed annulus to ^J{1 — pY + 77^ = 0.93(3), as shown in Fig.^J). 

The B^ — B^ lifetime difference is also an interesting quantity, but one for 
which theoretical progress is desirableM^ Lattice QCD can contribute by calcu- 
lating {Bs\(J)s)s-p(hs)s-p\Bs)- A preliminary result from JLQCD was reported 
some time ago: the ratio of the correct matrix element to the vacuu m-sa turation 
approximation was found to be Bs(TOb) = 0.86(3)(7) in the MS scheme.^^ 



5. Penguins and Associated Challenges 

5.1. Rare B Decays 

Lattice QCD could provide input to tests of the CKM model using radiative de- 
cays B -^ V"f. Belle has measured branching fractions for B -^ {p/tij)"fU^ Many 
new flavor models predict enhanced B -^ K*j decays. While LQCD input for the 
releva nt f orm factor is highly desirable, serious technical problems must first be 
solved.l^At the physical point, q'^ = 0, calculations done in the B rest frame have 
a final state hadron with momentum much greater than the inverse lattice spacing, 
ensuring lattice artifacts dominate. The challenge in working with small q^ with 
present methods, i.e. in the B rest frame, is that lattice artifacts dominate when 
the light meson's momentum becomes comparable to the inverse lattice spacing. 
Efforts are underway to formulate and employ lattice NRQCD in a frame where 
the B is moving, thus keeping the l ight meson's momentum small compared to the 
inverse lattice spacing! '^ ' ' * I ' v l 1 Vlff U jn addition, one again has to face questions 
regarding extrapolating lattice data through thresholds. 

Lattice QCD can play a more immediate role in rare B decays by calculating 
the form factors for B -^ K^*H'^(.~ decays. Branching fractions apparently agree 
with the Standard Model predictions,* ^ but these predictions rely on sum 

rule determinations for the form factors. As the experimental precision increases, 
the burden will turn to LQCD to reduce theoretical uncertainty. In fact, the LQCD 
calculation for the most important B — > Ki~^i~ matrix element (within the Stan- 
dard Model), {K\sjfj_b\B), is very similar to the B -^ ntv form factor calculations 
discussed above (Sec.|3|.^The only difference is that the mass of the strange quark 
should be held fixed as the mass of the spectator light quark is chirally extrapolated. 
B -^ K*£~^£~ would be worth exploring as well. 
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5.2. AB — Matrix Elements 

Lattice QCD can improve theoretical calculations of the B meson lifetimes through 
calculations of {B\Oab=o\B), where the O represents four AB = 4-quark 
operators?"^ Within the context of a certain flavor physics model, e.g. the CKM 
model, the experimental B meson (and Af,) lifetimes can be used to test the validity 
of the quark-hadron duality required to express the lifetimes in terms of those ma- 
trix elements. Two main challenges obstruct this calculation; these are familiar from 
LQCD calculations of (7r|OAS=i|^) for K -^ tttt. First, these operators mix with 
lower dimensional ones, and a power-law subtraction must be performed. Second, it 
is very difficult technically to calculate the contribution to the matrix element when 
the light quark and antiquark in 0{x) are contracted to form a propagator G{x, x). 
These challenges are not insurmo unt able and deserve to be met. An exploratory 
study was reported some time ago.l^ 

These calculations would have impact beyond the B meson lifetimes themselves. 
Again using quark-hadron duality and the optical theorem, LQCD calculations of 
the matrix elements of two AB ~ 4-quark operators would be helpful in improving 
the determination of \Vub\ from B -^ Xui~i'!^ 



6. Charmed Beauty 

Using the MILC configurations discussed above, a very precise prediction for the Be 
mass was recently obtained: m'g** = 6304 ± 12 +J^ MeV.EHlThis lattice QCD pre- 
diction has since been confirmed by experiment, m^^ = 6287.0 ± 4.8 ±1.1 MeV.!"^ 
The correlation function necessary for the calculation used a nonrelativistic action 
for the b quark and the Fermilab relativistic action for the c quark. Consequently, 
the agreement between LQCD and experiment further supports the use of both 
heavy quark formulations and improved staggered sea quarks. 

7. Conclusions 

The past 2 years have been especially fruitful for the interplay between lattice QCD 
and heavy flavor physics, with sharp improvements in LQCD and experimental 
results. Precise measurements of i? — > Trii/ decay combined with LQCD calculations 
of -B ^ TT form factors are making exclusive determinations of \Vub\ as precise 
as inclusive determinations, with prospects for further improvement. Employing 
an improved staggered quark action, LQCD calculations can be done with light 
quark masses inside the chiral regime. The resulting prediction of fs has since been 
confirmed with the first measurement of _B —> ri^ decay. Furthermore, removing 
chiral extrapolation and quenching errors from the ratio fss/fs combines with 
the first measurements of Arris to drastically improve the precision of \Vtd\- The 
orthogonality of the \Vtd\ and sin 2/3 constraints on (p, fj) creates a small target for 
other CKM constraints to hit. Most of us hope the Cabibbo-Kobayashi-Maskawa 
trio's aim will eventually fail, cracking the Standard Model. If not, the burden is on 
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the model builders to explain the remarkable craftsmanship of the CKM mechanism. 
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